PITHA 07/08 
MPP-2007-132 
UWThPh-2007-24 
PSI-PR-07-08 



00 
O 

o 
<^ 

u '. Radiative corrections to W-boson hadroproduction: 

! higher-order electroweak and supersymmetric effects 

. . SiLjA Brensing\ Stefan Dittmaier^'^, Michael Kramer\ Alexander Muck^'"^ 
^; 

^Institut fiir Theoretische Physik E, RWTH Aachen, 
D-52056 Aachen, Germany 

Max- Planck- Institut fiir Physik (Werner-Heisenberg-Institut), 
D-80805 Miinchen, Germany 

^ ' Faculty of Physics, University of Vienna, 

^ '. A- 1090 Vienna, Austria 

■ ^Paul Scherrer Institut, Wiirenlingen und Villigen, 

m 



o 



X 



Ch-5232 Villigen PSI, Switzerland 



Abstract: 

The high accuracy envisaged for future measurements of W-boson production at hadron 
cohiders has to be matched by precise theoretical predictions. We study the impact of electroweak 
radiative corrections on W-boson production cross sections and differential distributions at the 
Tevatron and at the LHC. In particular, we include photon-induced processes, which contribute 
at 0(0), and leading radiative corrections beyond 0{a) in the high-energy Sudakov regime 
and from multi-photon final-state radiation. We furthermore present the calculation of the 
complete supersymmetric next-to-leading-order electroweak and QCD corrections to W-boson 
hadroproduction within the MSSM. The supersymmetric corrections turn out to be negligible 
in the vicinity of the W resonance in general, reaching the percent level only at high lepton 
transverse momentum and for specific choices of the supersymmetric parameters. 



1 Introduction 



The Drell-Yan-like production of W bosons, 

PP/PP ^ W ^ li^iX, (1.1) 

is one of the cleanest hadron colhder processes with a large cross section at the Tevatron and at 
the LHC. Measurements near the W resonance allow for a precise determination of the W-boson 
mass and yield valuable information on the parton structure of the proton. Above resonance, 
the off-shell and high-energy tails of appropriate distributions are sensitive to the W width and 
offer the possibility to search for additional charged gauge bosons W. (See e.g. Refs. [1,2] and 
references therein.) 

The high experimental precision envisaged at the Tevatron Run II and specifically at the 
LHC has to be matched by precise theoretical predictions including QCD and electroweak ra- 
diative corrections. The QCD corrections arc completely known up to next-to-ncxt-to-leading 
order (NNLO) [3] and up to N^LO in the soft-plus- virtual approximation [4], with a remaining 
theoretical error for inclusive cross sections at the percent level or lower. The next-to-leading- 
order (NLO) QCD corrections have been matched with parton showers [5] and combined with a 
summation of soft gluon radiation (see e.g. Ref. [6]), which is particularly important to reliably 
predict the W transverse-momentum distribution. A theoretical study of the QCD uncertainties 
in the determination of the W cross section at hadron colliders has been presented in Ref. [7]. 

In this article we mainly focus on the electroweak corrections. Given, for example, the 
anticipated experimental accuracy in the W-boson mass measurement of 10— 20MeV at the 
LHC, the inclusion of electroweak corrections beyond final-state radiation is mandatory as their 
omission would induce a systematic error in the Mw determination of 0(10 MeV) (see e.g. 
Ref. [1]). The complete NLO 0{a) corrections to the parton process gg' — ^ W — > lui have 
been calculated by several groups [8-11]. A tuned comparison of cross sections and differential 
distributions has shown good agreement between the various calculations [1,12]. Photon-induced 
processes 'jq l^iq' , which also contribute at 0{a), have been computed recently in Refs. [13, 14] . 
The particular importance of final-state photon radiation for the W-boson mass determination 
demands a treatment that goes beyond 0{a). Such multi-photon effects have been studied 
in Refs. [15,16] and matched to the NLO 0{a) calculation in Ref. [11]. First steps towards 
combining QCD and electroweak higher-order effects have been taken in Ref. [17]. We note that 
the 0{q) corrections to hadronic production of on-shell W bosons at large transverse momenta, 
PP/pP ^ W -I- jet, have recently been presented in Ref. [18] for the Standard Model (SM) 
and in Ref. [19] in the minimal supersymmetric Standard Model (MSSM). For related work on 
electroweak corrections to Z-boson hadroproduction we refer to Refs. [1, 12] and further references 
therein. 

In this paper we improve on our previous 0{a) calculation [9] for pp/pp -^W ^ luiX by in- 
cluding photon-induced processes and multi-photon final-state radiation in the structure- function 
approach. We also discuss the impact of leading electroweak effects beyond 0{a), specifically 
the Sudakov logarithms that arise in the high-energy regime, and present the calculation of the 
complete supersymmetric 0{a) electroweak and 0{as) strong corrections within the MSSM. The 
purpose of the MSSM calculation is to establish that the impact of virtual supersymmetric par- 
ticles on the cross section prediction is small and does not spoil the status of single- W-boson 
production as one of the cleanest SM candles at hadron colliders. 

The paper is organized as follows. In Section 2.1 we briefly summarize the calculation of 
the complete 0{a) corrections as originally presented in Ref. [9]. In the following sections, 
we describe the inclusion of leading electroweak effects beyond 0{a), the calculation of the 
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photon-induced processes, and the summation of multiple emission of collinear photons off the 
final-state lepton. Numerical results for W-boson production at the Tevatron and the LHC 
including the complete 0{a) and 0(as) NLO corrections and the higher-order Sudakov and 
photon radiation effects are presented in Section 2.6. There, we also compare our results on 
multi-photon final-state radiation in the structure-function approach with results based on the 
parton-shower approach [1]. In Section 3 we describe the calculation of the MSSM corrections 
and discuss their numerical impact on the W cross section and distributions. We conclude in 
Section 4. Finally, the Appendix provides details on the scenarios of the supersymmetric models 
under consideration. 



2 Higher-order electroweak effects 

In the next section, we briefly review the calculation of the complete 0{a) corrections to the 
parton process qq' — > W — > Iz/; presented in Ref. [9] which is now augmented by an extension 
of the dipole subtraction method [20] that allows us to calculate non-collincar-safc obscrvables 
also within the subtraction approach. In the subsequent sections, we then discuss the leading 
electroweak effects and the choice of couplings before we present those parts of the calculation 
that are new and that have not been discussed in Ref. [9] . We present our results for the specific 
case of W"*" production, pp/pp — W+ — > Z+f/X. 

2.1 Electroweak 0(cy.) corrections 

We consider the parton process 

n{pu) + d{pd) ^ Mkn) + l'^{ki) [+l{k)], (2.1) 

where u and d generically denote the light up- and down-type quarks, ti = u, c and cZ = d, s. 
The lepton I represents I = e,fi. The momenta of the particles are given in brackets, and the 
Mandelstam variables are given by 

s = {Pu+Pd?, i={pd-kif, u={pu-hf- (2.2) 

We neglect the fermion masses m^, m^, mi whenever possible, i.e. we keep these masses only 
as regulators in the logarithmic mass singularities originating from collinear photon emission or 
exchange. In lowest order the scattering amplitude reads 

= f''^^'"-""^ .--M^ + iMwrw(^) t^'^'^'^"-^'^ ' ^'-'^ 

with an obvious notation for the Dirac spinors Vd, etc., and the left-handed chirality projector 
oj- = ^(1 — 75). The electric unit charge is denoted by e, the weak mixing angle is fixed by 
the ratio = 1 — = M-^/Af| of the W- and Z-boson masses Mw and Mz, and Vud is the 
CKM matrix element for the ud transition. The different choices for the fine-structure constant 
a = e^/(47r) in the squared matrix element are discussed below. 

Strictly speaking, Eq. (2.3) already goes beyond lowest order, since the W-boson width rw(s) 
results from the Dyson summation of W self-energy insertions. Renormalizing the W mass and 
width in the on-shell scheme, the Dyson summation directly leads to a running width, i.e. 

rw(^)L,. = rw^. (2.4) 
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On the other hand, defining W mass and width from the location of the pole of the W propagator 
(with momentum transfer p) in the complex plane, naturally leads to a constant width, i.e. 

rw(5)L_t = rw. (2.5) 

In this work, we employ the fixed- width approach (2.5). This has to be kept in mind when using 
the results in precision determinations of the W-boson mass from the W resonance. 

The virtual one-loop corrections comprise contributions of the transverse part of the W self- 
energy Y^y^ , corrections to the two Wdu and Wvil vertices, box diagrams, and counterterms. 
Our calculation of these corrections is described in Ref. [9] in detail. In particular, the complete 
expressions for the vertex and box corrections are provided in the Appendix of that reference. For 
this paper we recalculated the one-loop effects employing the packages FeynArts [21], FormCalc, 
and LoopTools [22], supplemented by the loop integrals with singularities on the W resonance, 
which are discussed below. In both calculations, ultraviolet divergences are treated in dimen- 
sional regularization, and the infrared (IR) singularity is regularized by an infinitesimal photon 
mass m-y. The actual calculation is performed in 't Hooft-Feynman gauge using on-shell renor- 
malization. 

The complete one-loop amplitude M.i can be expressed in terms of a correction factor 5^'^ 
times the lowest-order matrix element, 

Ml = (5™*7Wo. (2.6) 

Thus, in 0{a) the squared matrix element reads 

|7Wo + VWiP = (l + 2Re{5"'^*})|7WoP + -.., (2.7) 

so that the Breit-Wigner factors are completely contained in the lowest-order factor |AloP- Note 
that the Dyson-summed imaginary part of the W self-energy, which appears as T^f^{s) in A^O; 
is not double-counted, since only the real part of enters Re{(5'^^'^*} in 0{a). Despite the 
separation of the resonance pole (s — M^)~^ from the correction factor 5^^"^ in (2.6), 5™^ still 
contains logarithms ln(s — M-^+ie) that are singular on resonance. Since these singularities would 
be cured by a Dyson summation of the W self-energy inside the loop diagrams, we substitute 

ln(s-M^ + ie) ^ ln(s - + iMwTw) (2.8) 

with a fixed width everywhere. The substitution (2.8) does not disturb the gauge-invariance 

properties of the one-loop amplitude ^Al. 

In Ref. [9] (Section 2.2) three different input-parameter schemes have been specified for the 
choice of the electromagnetic coupling constant a that is taken as SM input together with the 
particle masses. According to the choice of a, the schemes are called "a(0)", "a;(M|)", and 
"G^" schemes, where a is set to ac^ = \/2G^M^(l — M^/M^)/tt in the scheme. In the 
a(0) scheme, the charge renormalization constant 6Z(, in the counterterm <5^yjj/ of the Wff 
vertex contains logarithms of the light-fermion masses which are related to the running of the 
electromagnetic coupling a{Q'^) from Q = to a high-energy scale. In the a(M|) scheme the 
counterterm changes to 

Swff'liMi) = ^wff'Uo) - ^Aa(M|), (2.9) 

where 

Aa{Q') = Uff,{0) - Re{Uff,{Q')}, (2.10) 
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with ^f^i denoting the photonic vacuum polarization induced by all fermions other than the 

top quark. In contrast to the a(0) scheme the counterterm ^wff'\a{M'^) ^o^^ involve light 
quark masses, i.e. the large logarithmic corrections are absorbed in the lowest order by using the 
appropriate numerical value for q;(M|). In the scheme, the transition from a(0) to is 
ruled by the quantity Ar which is deduced from muon decay, 

= ^^M^w4 ^ ^(o)(i + + 0(^3)_ (2.11) 

TT 

Therefore, the counterterm S^jrj^, reads 

Since Aa(Af|) is implicitly contained in Ar, the large fermion-mass logarithms are also absorbed 
in the lowest order in the scheme. The different input-parameter schemes are further discussed 
in Section 2.2. 

The calculation of the real-photonic corrections is described in Ref. [9] in detail, both for a 
running and a fixed W-boson width. In order to respect electromagnetic gauge invariancc, the 
coupling of the photon to W bosons has to be adapted when a running width is used. In the 
fixed- width approach, as used in this paper, no modification of couplings is needed. 

The helicity amplitudes for the radiative process are explicitly given in Section 2.4 of Ref. [9]. 
The phase-space integral over the real-photonic matrix elements diverges in the soft and collinear 
regions logarithmically if the photon and fcrmion masses are set to zero. To properly combine the 
soft and collinear singularities with the corresponding virtual corrections three different methods 
are applied: two variants of phase-space slicing and the dipole subtraction method. 

The dipole subtraction approach as formulated in Ref. [23] can only be applied for collinear- 
safe observables, i.e. observables for which selection cuts are blind to the distribution of mo- 
menta in collinear lepton— photon configurations. This can be achieved by "photon recombina- 
tion", where leptons and sufficiently collinear photons are treated as one quasi-particle (see also 
Section 2.6.2). For these observables the KLN theorem [24] guarantees that logarithms of the 
fermion masses are absent in the corrections. For muons in the final state it is, however, experi- 
mentally possible to separate collinear photons from the lepton, i.e. to observe so-called "bare" 
muons. Hence, the resulting cross sections are not collinear safe and the corresponding collinear 
singularities show up as logarithms of the small lepton masses. 

In Ref. [9], only the slicing variants were able to deal with bare leptons, while the application 
of the subtraction approach was still restricted to collinear-safe observables. In this work, we 
employ an extension [20] of the subtraction formalism which allows one to calculate cross sections 
for bare leptons, i.e. cross sections defined without any photon recombination. The respective 
results of the slicing and subtraction methods are in good numerical agreement both for the 
inclusive and the bare muon case. 



2.2 Leading electroweELk effects and choice of couplings 

As described in Ref. [9], the relative corrections in the various input-parameter schemes differ 
by constant contributions proportional to Aa(M|) ~ 6% and Ar ~ 3%, which quantify the 
running of the electromagnetic coupling from = to = M| for q;(M|) and the radiative 
corrections to muon decay for the scheme, respectively. The bulk of Ar is contained in 
Aq;(M|) — c^A/j/s^, where Ap ^ 1% is the universal correction oc G^ml to the p parameter. 
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The Gfj_ scheme is distinguished from the two other schemes because the corrections to 
charged-current four-fermion processes do not contain large contributions from Aa(M|) or Ap 
anymore, i.e. these universal renormalization effects are completely absorbed into the leading- 
order (LO) amplitude. For the 0{a) corrections this property has already been pointed out in 
Ref. [9], but it also holds at 0{a'^). More precisely, in the G^j scheme no contributions propor- 
tional to Aa(Af|)" (for any positive integer n), proportional to Aa(Af|)Ap, and proportional to 
Ap^ appear. Based on the arguments given in Ref. [25] this feature was explicitly worked out in 
Ref. [26] (Section 3) for the related process of charged-current neutrino deep-inelastic scattering; 
these results obviously apply also to charged-current Drcll-Yan scattering via crossing symmetry. 

Before concluding that the theoretical uncertainties in the scheme from missing corrections 
beyond 0{a) are smaller than in the two other schemes, we shall inspect other known universal 
dominant corrections at 0{a). Besides the renormalization effects discussed above, the dominant 
0{a) corrections, up to moderate parton scattering energies, are due to final-state radiation off 
the charged lepton, at least for bare leptons where the enhancement by the large mass logarithm 
oc a In mi is present. The inclusion of these contributions beyond 0{a), which are due to collinear 
multi-photon emission, is described in Section 2.5 below. Here we merely point out that the ap- 
propriate coupling constant a entering the relative correction is a{0), because it accounts for the 
emission of photons with = 0. Thus, when adopting an input-parameter scheme other than 
the a(0) scheme, one should nevertheless use a(0) to multiply the dominating universal lepton- 
mass logarithms. Our specific implementation of this procedure is described in Section 2.5. 

At high parton scattering energies and high transverse lepton momenta, electroweak correc- 
tions are dominated by soft and/or collinear gauge-boson exchange. The soft effects, known as 
Sudakov logarithms, induce powers of Q;ln^(s/M^), with subleading soft /collinear contributions 
involving lower powers in the logarithm. In 0{a), these enhanced electroweak effects drive the 
relative corrections to ^ —30% at lepton transverse momenta of about ITeV at the LHC [9]. 
It is therefore desirable to control this kind of corrections beyond 0{a). We elaborate more on 
this issue in the next section. Here we just point out that among the considered input-parameter 
schemes the scheme should be appropriate to fix a for the leading high-energy logarithms, 
which are of weak origin. 

Following the above arguments, we employ the G^ scheme in this work, modified only by the 
change of ao^ to a(0) in the leading part of final-state radiation. This procedure is expected to 
be most robust with respect to further corrections beyond 0{a). 

2.3 Leading weeik corrections in the Sudaikov regime 

For single-W production at large lepton transverse momenta or W transverse masses, the parton 

kinematics is restricted to the Sudakov regime, characterized by large Mandelstam parameters s, 
\u\ 3> -M^. The structure of electroweak corrections beyond 0{a) in this high-energy regime 
has been investigated in some detail by several groups in recent years (see e.g. Refs. [27-34] and 
references therein). 

As described for example in Refs. [32,34], the leading electroweak logarithmic corrections, 
which are enhanced by large factors L = ln(s/M^), can be divided into an SU(2)xU(l)- 
symmetric part, an electromagnetic part, and a subleading part induced by the mass differ- 
ence between W and Z bosons. The last part does not contribute to corrections oc (ccL^)" and 
is neglected in the following. The leading (Sudakov) logarithms oc (aL^)" of electromagnetic 
origin cancel between virtual and real (soft) bremsstrahlung corrections; for the subleading log- 
arithms such cancellations should strongly depend on the observable under consideration. The 
only source of leading logarithms is, thus, the symmetric electroweak (sew) part, which can be 
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characterized by comprising W bosons, Z bosons, and photons of a common mass Mw- Using 
this mass assignment, the one- loop correction to the squared amphtude can be obtained 

rvirt 



by expanding the full result for the virtual correction S™^ (given in Appendix A of Ref. [9]) for 

f2 



large s, \u\ M^. The exphcit result can be written as 



dig = ^ {-L'C^^c + LCr},c} (2.13) 



with factors 



In — + In 



-u 



+ 2ly Y 



-^l.CC - ^ ^ 73 ' "^l,CC - --2- 

(2.14) 

which have been introduced in Section 8.4.2 of Ref. [34]. Here = 1/3 and Yi,^^ = —1 are 
the weak hypercharges of the corresponding left-handed particles. In Eq. (2.13) wc did not only 
include the leading Sudakov logarithms cx aL^, but also the related "angular-dependent" con- 
tributions oc aL\n{—t/s) or aLhi{—u/s). Our explicit 0{a) result is in agreement with the 
general results presented in Refs. [32,34], where the corresponding corrections are also given at 
the two- loop level. These 0{a'^) corrections can be obtained from the 0{a) result by an appro- 
priate exponentiation [30]. For the leading "sew" corrections (including a^L^, ct^L^ ln(— t/s), 
and Li^ \n.{—u / s) terms) this exponentiation simply reads [34] 

\M\^ ~ |Mo|' exp {<5W } = \Mo\^ (l + 5^1 + 4i + . . . ) (2.15) 

with 

4i ={^)\ \^'(CfZc? - L^Cf^cCrhc } • (2-16) 

However, in the case of neutral-current-induced fermion-antifermion scattering processes it 
was observed [33] that large cancellations take place between leading and subleading logarithms. 
In view of this uncertainty, we do not include the two-loop high-energy logarithms in our full pre- 

(2) 

dictions. Instead, we evaluate the leading two-loop part ^sew 

as a measure of missing electroweak 
corrections beyond 0{a) in the high-energy Sudakov regime. 

Moreover, since the electroweak high-energy logarithmic corrections are associated with vir- 
tual soft and/or collinear weak-boson or photon exchange, they all have counterparts in real 
weak-boson or photon emission processes which can partially cancel (but not completely, see 
Ref. [28]) the large negative corrections. To which extent the cancellation occurs depends on 
the experimental possibilities to separate final states with or without weak bosons or photons. 
This issue is discussed for example in Refs. [35, 36] . The numerical analysis presented in Ref. [36] 
demonstrates the effect of real weak-boson emission in the distributions in the transverse lepton 
momentum px,z and in the transverse mass M^^^^^i of the W boson (as e.g. defined in Section 2.6.3 
below). For W"*" production at the LHC, at M^^^^i = 2TeV the electroweak corrections are re- 
duced from about —26% to —23% by weak-boson emission. At p^^i = 1 TeV the corresponding 
reduction from about —28% to —17% is much larger, however, the bulk of these emission effects 
is not due to soft /collinear weak-boson emission, but due to recoil effects in real W^W^ and 
W+Z production. This explicitly illustrates the sensitivity of weak-boson emission effects to the 
details of experimental event selection, in particular, how single-W production is separated from 
di-boson production. 
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2.4 Photon-induced processes 

The 0{a) corrections to the parton cross section qq' W+ I'^vi contain coUinear singularities 
from photon radiation off the initial-state quarks which arc absorbed by mass factorization [9-11]. 
For a complete and theoretically consistent analysis, the absorption of the 0{a) coUinear singu- 
larities into quark distributions has to be complemented by both the inclusion of 0{a) corrections 
to the parton distribution functions (PDFs) and by the inclusion of the 0{a) partonic subpro- 
cesses 7U — > I'^uid and 7d — > I'^vin, which were first calculated in Ref. [13]. At the time of 
our previous study [9] a complete analysis of 0{q) corrections to PDFs, which involves 0{a) 
corrections to the DGLAP evolution and to the fit of experimental data, was not available. 
Meanwhile, the MRSTQED2004 [37] PDF parametrization, however, includes these 0{a) cor- 
rections and provides a photon density within protons and antiprotons. It is thus now possible 
to perform a complete 0{a) analysis and to properly include the photon- induced subprocesses 
in the hadronic cross-section prediction. 

The coUinear photon splitting into two massless quarks in the subprocesses 7U — I'^vid and 
7d ^ I^uiVl also leads to a mass singularity. As usual, this divergence is removed by mass 
factorization. Including both the divergences from collinear photon splitting and the divergence 
due to collinear photon emission from initial-state quarks in qq' l^i^n, mass factorization 
implies a redefinition of the quark densities according to [26] 



dz 
z 



x<!ln(4 



, (X \ r^2 ^ 
Pqq{z) 



- Pqq(z)(21n(l-z) + l)J^ + Cqq(z) 
£|;,g,,,)3gjii{l„(|)p„(.)+C„(.)}, 



(2.17) 

where Qq is the electric quark charge, is the small quark mass used as a regulator, and y^p 
denotes the QED factorization scale which is identified with the QCD factorization scale. The 
factor 3 in the second line stems from the splitting of the photon into qq pairs of different color. 
Furthermore, 

W = ^, P,,{z)=z^^(\-zf (2.18) 



1 



are the quark and photon splitting functions, respectively, and Cgg, C^-y the coefficient functions 
specifying the factorization scheme. Following standard QCD terminology one distinguishes MS 
and DIS schemes defined by 



Cf^{z) = Cf^{z) = 0, 



Pqq{z) In 



\-z 



3\ 9-I-5Z 
4/ ^ 4 



(2.19) 



J + 



CDIS(^) = P^^(^)ln 



1 



- -h 8z - 1 . 



In our numerical analysis we employ the MRSTQED2004 parton distribution functions. Note 
that photon radiation off incoming quarks was ignored in the F2 fit to HERA data in the 
MRSTQED2004 PDF determination. Therefore, the MRSTQED2004 PDFs are defined in the 
DIS scheme for the factorization of QED effects, i.e. not in the MS scheme as frequently done in 
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the past (see also Ref. [26]). For the factorization of QCD effects the PDFs are defined in the 
MS scheme as usual. 

To extract the colhnear divergence from the squared matrix element for the photon-induced 
processes we use an extension [20] of the dipole subtraction technique, which has been formulated 
to treat the collinear splitting of photons into light fermions 7 ff. 

2.5 Multi-photon final-state radiation 

The emission of photons collinear to the outgoing charged lepton leads to corrections that are 
enhanced by large logarithms of the form aln{mf /Q'^) with Q denoting a characteristic scale of 
the process. The KLN theorem [24] guarantees that these logarithms cancel if photons collinear 
to the lepton are treated fully inclusively. However, since we apply a phase-space cut on the 
momentum of the outgoing lepton, contributions enhanced by these logarithms survive if the 
momentum of the bare lepton is considered, i.e. if no photon recombination is performed. While 
the concept of a bare lepton is not realistic for electrons, it is phenomenologically relevant for 
muon final states. 

The first-order logarithm aln{mf /Q"^) is, of course, contained in the full 0{a) correction, so 
that Q is unambiguously fixed in this order. However, it is desirable to control the logarithmically 
enhanced corrections beyond 0{a). This can be done in the so-called structure-function approach 
[38], where these logarithms are derived from the universal factorization of the related mass 
singularity. The incorporation of the mass-singular logarithms takes the form of a convolution 
integral over the LO cross section ctq. 



where the step function Ocut is equal to 1 if the event passes the cut on the rescaled lepton 
momentum zki and otherwise. The variable z is the momentum fraction describing the lepton 
energy loss by collinear photon emission. Note that in contrast to the parton-shower approaches 
to photon radiation (see e.g. Refs. [11, 15, 16]), the structure-function approach neglects the pho- 
ton momenta transverse to the lepton momentum. 

For the structure function Tjf"{z, Q^) we take into account terms up to 0{a^) improved by 
the well-known exponentiation of the soft-photonic parts [38], 




(2.20) 




(2.21) 



8 



with and T{y) denoting Eulcr's constant and the Gamma function, respectively. The large 
logarithm is contained in the variable 



2a(0) 



TT 



1 



(2.22) 



Here, a(0) is the fine-structure constant defined in the Thomson limit. The parts solely pro- 
portional to a power of (3i correspond to collinear (multi-)photon emission off the lepton, the 
exponential factor describes resummed soft-photonic effects. The non-logarithmic term "—1" in 
f3i accounts for a non-singular universal soft-photonic correction. 

Technically, we add the cross section (2.20) to the one-loop result and subtract the LO and 
one-loop contributions 



^LLiFSR= / dcro(pu,Pd;^i/nfcO dz ^d{l - z) + r\f^''^ (z, Q'"^)^ @cut{zki), 



(2.23) 



contained in (2.20) in order to avoid double counting. The one-loop contribution to the structure 
function reads 



4 V 1-z 



(2.24) 



More precisely, we adapt the value of a in T^^'^{z, Q^) to the chosen input parameter scheme. 
Thereby, we introduce an additional higher-order contribution (a(0) — a) In^mf/Q'^) so that the 
aln{mf /Q^) contribution to the full 0{a) correction is subtracted exactly. Hence, the procedure 
of adding higher-order final-state radiation changes also the value of a in the Q;ln(m^/Q^) term 
to a(0) which is the appropriate coupling for real-photonic effects. 

The uncertainty that is connected with the choice of enters in ©(a^), since all 0{a) 
corrections, including constant terms, are taken into account. As default we choose the value 



with ^ 



1. In order to quantify the scale uncertainty, we vary ^ between 1/3 and 3. 



(2.25) 



2.6 Numerical results 

2.6.1 Input parameters and setup 

The relevant SM input parameters are 





= 1.16637 X 10-5 GeV-^ 


a{0) 


= 1/137.03599911, 


as(Mz) 


= 0.1189, 


Mw 


= 80.403 GeV, 


Tw 


= 2.141 GeV, 






Mz 


= 91.1876 GeV, 


Mh 


= 115 GeV, 






me 


= 0.51099892 MeV, 


rUf, 


= 105.658369 MeV, 


rur 


= 1.77699 GeV, 


rriu 


= 66 MeV, 


rUc 


= 1.2 GeV, 


nit 


= 174.2 GeV, 


rud 


= 66 MeV, 


nis 


= 150 MeV, 


mb 


= 4.6 GeV, 


iKdl 


= 0.974, 


iKsl 


= 0.227, 






IKdl 


= 0.227, 




= 0.974, 







(2.26) 
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which essentially follow Ref. [39]. The masses of the light quarks are adjusted to reproduce 
the hadronic contribution to the photonic vacuum polarization of Ref. [40]. The CKM matrix 
is included via global factors in the partonic cross sections for the different initial-state quark 
flavours. Within loops the CKM matrix is set to unity. 

As explained in Section 2.2, we adopt the scheme (up to the modification of a in the final- 
state radiation) , where the electromagnetic coupling a is set to . The charge renormalization 
constant, which contains logarithms of the fermion masses, drops out in the scheme so that 
our results are practically independent of the light-quark masses. We keep finite light-quark 
masses in closed fermion loops, their nTimcrical impact is however extremely small. The W- 
boson resonance is treated with a fixed width without any running effects. 

The 0(a)-improved MRSTQED2004 set of PDFs [37] is used throughout. The QCD and 
QED factorization scales are identified and set to the W-boson mass Mw- 

2.6.2 Phase-space cuts and event selection 

For the experimental identification of the process pp/pp W+ l^i^iX we impose the set of 
phase-space cuts 

PT,/>25GeV, |^T>25GeV, |??;| < 2.5, (2.27) 

where pt,z and rn arc the transverse momentum and the rapidity of the charged lepton l'^ , re- 
spectively, and i^T = P'T,ui is the missing transverse momentum carried away by the neutrino. 
Note that compared to our previous study [9] we have changed the rji-cuX, to \r]i\ < 2.5, which is 
a more realistic estimate of the experimental charged-lepton coverage at the LHC. The identifi- 
cation cuts are not collinear safe with respect to the lepton momentum, so that observables in 
general receive corrections that involve large lepton-mass logarithms of the form Q;ln(m;/Mw). 
This is due to the fact that photons within a small collinear cone around the charged-lepton 
momentum are not treated inclusively, i.e. the cuts assume a perfect isolation of photons from 
the charged lepton. While this is (more or less) achievable for muon final states, it is not realistic 
for electrons. In order to be closer to the experimental situation for electrons, the following 
photon recombination procedure is applied: 

1. Photons with a rapidity |ry^| > 3, which are close to the beams, are considered part of the 
proton remnant and are not recombined with the lepton.^ 

2. If the photon survived the first step, and if the resolution Ri^ = ^/^(r?i^-^?TO^~+~0f^ is smaller 
than 0.1 (with (f)ij denoting the angle between lepton and photon in the transverse plane), 
then the photon is recombined with the charged lepton, i.e. the momenta of the photon 
and of the lepton I are added and associated with the momentum of I, and the photon is 
discarded. 

3. Finally, all events are discarded in which the resulting momentum of the charged lepton 
does not pass the cuts given in (2.27). 

While the electroweak corrections differ for final-state electrons and muons without photon re- 
combination, the corrections become universal in the presence of photon recombination, since 
the lepton-mass logarithms cancel in this case, in accordance with the KLN theorem. Numerical 
results are presented for photon recombination and for bare muons. 

^Note that collinear safety requires that the jr;^| cut must be larger than the lepton identification cut on \rii\ 
to avoid events where an almost collinear lepton-photon pair is not recombined because \rji\ < 2.5 and \rij\ ^ 2.5. 
It turns out, however, that the numerical diiTerence between choosing \ri~i\ > 3 and \ri^\ > 2.5 is negligible. 
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pp l+viX at = 14TeV 



PT,i/GeV 


25-00 


50-00 


lOO-oo 


200-00 


500-00 


1000-00 


cTo/pb 


4495.7(2) 


27.589(2) 


1.7906(1) 


0.18128(1) 


0.0065222(4) 0.00027322(1) 




-2.9(1) 


-5.1(1) 


-8.6(1) 


-13.2(1) 


-23.4(1) 


-34.7(1) 




-1.8(1) 


-2.6(1) 


-6.1(1) 


-10.3(1) 


-19.5(1) 


—29.5(1) 




0.065(1) 


4.7(1) 


12.3(1) 


17.1(1) 


16.7(1) 


13.5(1) 


e(2) /q/ 
"Sudakov//^ 


—0.0002 


—0.023 


—0.082 


0.057 


1.3 


3.8 


(>multi-7/ /O 


„ , n + O 03 


„ +0 08 


„ r,r-,+0 06 
0.27;^:^° 


^ „T +0 06 

0-31Io:o6 


n 41 +0-08 

^•^■'--0.07 


-0.09 




-2.7(1) 


0.0(1) 


4.0(1) 


4.3(1) 


-6.3(1) 


-20.6(1) 




-1.7(1) 


2.1(1) 


6.2(1) 


6.9(1) 


-2.7(1) 


-16.0(1) 




-2.7(1) 


812(1) 


784(1) 


814(1) 


611(1) 


399(1) 




-2.8(1) 


793(1) 


685(1) 


607(1) 


323(1) 


127(1) 



Table 1: Integrated LO cross sections ao for W+ production at the LHC for different 
ranges in px,/ and corresponding relative corrections S in the SM. 



2.6.3 Cross sections and distributions for pp — > W+ — > I'^uiX at the LHC 

We first consider W"*" production at the LHC, i.e. a pp initial state with a centre-of-mass (CM) 
energy of ^/s = 14 TeV. 

In Tables 1 and 2 we present the LO cross section gq and various types of electroweak 
corrections 5, defined relative to the LO cross section by a = x [1+5). The results are 
shown for different ranges in the transverse momentum of the charged Icpton, Pt,«) and in the 
transverse mass of the two final-state leptons, M^^^vii = y^2pT,i|^T(l — cos(/)^,,i), where (t)vii is the 
angle between the lepton and the missing momentum in the transverse plane. 

For reference, we first update the 0{a) NLO corrections for W-boson hadroproduction 
through the parton process qq' l^ui [9]. The results are given for bare muon final states 

{Sgg '^'^) and with photon recombination applied {5^'^^). As explained above, the mass-singular 

corrections oc aln{mn/My^) present in Sgg '^'^ cancel if the photon is recombined, rendering the 
corresponding correction smaller. At large p^^i and M^^^^i the electroweak corrections are 
dominated by the 0{a) Sudakov logarithms discussed in Section 2.3. Note that the relative 0{a) 
corrections 5qq presented in Tables 1 and 2 are not very sensitive to the choice of the PDF and 
the choice of the rji cut and thus agree very well with our previous numerical results presented 
in Ref. [9]. 

The 0{a) corrections originating from the photon-induced processes (Section 2.4) are not 
included in 5qq, but are shown separately as J-y^ in Tables 1 and 2. They are enhanced at large 
Pt,/ because of a new type of contribution where the incoming photon couples to a W boson 
that is exchanged in the t-channel. The photon-induced processes could in principle be used to 
extract information on the photon content of the proton. However, they are overwhelmed by QCD 
corrections and QCD uncertainties which strongly affect the pT,i spectrum (see the discussion 
of NLO QCD corrections below). If, on the other hand, one considers the distribution in the 



11 



pp l+viX at = 14TeV 



MT,.,i/GeV 


50-00 


100-00 


200-00 


500-00 


lOOO-oo 


2000-00 


(7o/pb 


4495.7(2) 


27.589(2) 


1.7906(1) 


0.084697(4) 0.0065222(4) 0.00027322(1) 




-2.9(1) 


-5.2(1) 


-8.1(1) 


-14.8(1) 


-22.6(1) 


-33.2(1) 




-1.8(1) 


-3.5(1) 


-6.5(1) 


-12.7(1) 


-20.0(1) 


—29.6(1) 




0.052(1) 


0.12(1) 


0.25(1) 


0.37(1) 


0.39(1) 


0.36(1) 


c(2) yq/ 

^SudaJiov/ 


—0.0002 


—0.023 


—0.082 


0.21 


1.3 


3.8 


Omulti-7/ /o 


„ _ n + O 03 


o.2oi^:^^ 


„ ^ „-i_n 03 


"-'•-'-^-0.03 




^•"^4-0.05 




-2.7(1) 


-4.9(1) 


-7.7(1) 


-14.2(1) 


-22.0(1) 


-32.5(1) 




-1.7(1) 


-3.4(1) 


-6.3(1) 


-12.3(1) 


-19.6(1) 


-29.3(1) 




-4.2(1) 


23.2(1) 


26.6(1) 


19.1(1) 


4.7(1) 


-18.5(1) 


<;H=Mt,W /O/ 

'^QCD / /o 


-4.4(1) 


22.5(1) 


24.0(1) 


12.6(1) 


-6.2(1) 


-34.6(1) 



Table 2: Integrated LO cross sections (Tq for W+ production at the LHC for different 
ranges in M^^i^^i and corresponding relative corrections 6 in the SM. 



transverse mass M^^^^^i, which is much less sensitive to QCD effects, the impact of 5^q is below the 
percent level. The results for S^g in Tables 1 and 2 are in good agreement with those presented 
in Refs. [13, 14] . Note that in Ref. [14] the mass factorization of the collinear j ^ qq splitting is 
performed in the MS scheme, while the photon distribution of the MRSTQED2004 PDF set is 
defined in the DIS scheme as argued above. It turns out, however, that the difference between the 
calculations in the MS and in the DIS scheme is numerically negligible for hadronic cross sections. 
We find that the 0{a'^) high -energy Sudakov logarithms calculated in Section 2.3, labeled 

(2) 

'^Sudakov Tables 1 and 2, have a small impact on the cross-section prediction, below 5% even 
for a transverse lepton momentum p^^i in the TeV range. 

The corrections due to multi-photon final-state radiation beyond 0{a) (see Section 2.5) are 
shown as 5muiti-7 in the tables. Only the genuine higher-order photon effects are included in 
'^muiti-77 the one-loop contribution is subtracted. We show (5inuiti-7 for the central scale choice 
Q = y/I with an uncertainty estimate obtained from varying the scale Q between Q = sVI (upper 
number) and Q = V§/3 (lower number). Multi-photon final-state radiation beyond 0{a) has a 
very small effect on the cross sections displayed in Tables 1 and 2. The largest part of this small 
contribution is in fact due to the change of the coupling constant a from ac"^ to a(0) in the relative 
0{a) correction (see also Section 2.6.5). However, as we shall discuss below, the contribution 
from multi-photon final-state radiation reaches the percent level near the W resonance and has 
thus a significant impact on a precision determination of the W mass. A more detailed analysis 
is needed to quantify the corresponding shift in the determination of Mw (cf. Refs. [1,15,41]). 

For convenience we combine the above results and display our best estimate for i) the elec- 

troweak corrections for muon final states (^ew*^) which includes the 0{a) correction to the qq' 
initial states, the corrections due to the photon-induced processes, and the multi-photon final- 
state radiation corrections with the scale choice Q = s/I, and ii) the total electroweak correction 
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pp l+viX at = 14TeV 



PT,//GeV 


25-00 


50-00 


lOO-oo 


200-00 


500-00 


lOOO-oo 


'^(77, veto /% 


0.025(1) 


1.2(1) 


0.049(1) 


0.043(1) 


0.042(1) 


0.042(1) 


jrtJ.=Mw IW 

"QCD.veto/ 


-7.3(1) 


454(1) 


6.4(1) 


-15.8(1) 


-51.9(1) 


-85.0(1) 



Table 3: Relative corrections for production at the LHC from the photon- induced 
processes and from the NLO QCD calculation with a jet veto imposed. We require the 
additional parton to be produced at PT,jet < 50 GeV. 



for final states with photon recombination (^g^, which combines the Oia) corrections from grg' 
and 7g initial states and which is not sensitive to multi-photon final-state radiation. Because 
of the theoretical uncertainty in evaluating the higher-order weak corrections in the high-energy 
regime (sec the discussion in Section 2.3), we do not include the leading two-loop Sudakov loga- 

(2) 

rithms (^gydakov ^^'^ '^^^^ estimate of the electroweak corrections. 

For comparison, we have also calculated the NLO QCD corrections, evaluated with two 
different choices for the renormalization and factorization scales, fiR = fip = -^w ('^qcd^^) ^^'^ 

HR = HF = Mt,w with M| w = M^+p^^^ (Sqcb)^ where PT,w is the transverse momentum 
of the W boson to be evaluated on an event-by-event basis. We have compared our NLO QCD 
results for Sq^^^ with those obtained from MCFM [42] and find good agreement. In the QCD 
case, there is of course only initial-state radiation and thus no technical need for a recombination 
procedure of the charged lepton with a possible additional jet. Rather, one should employ a 
separation cut between lepton and jet to allow for a clean event selection. However, since we 
only want to give a rough estimate of the size of QCD effects for comparison with the electroweak 
corrections, for simplicity, we do not include any sort of separation cut. As indicated above, the 
QCD corrections are extremely large at large px,; so that the electroweak corrections to the px,; 
distribution are overwhelmed by QCD uncertainties. The M^^^y^i distribution, on the other hand, 
is much less sensitive to higher-order QCD effects. It is invariant under transverse boosts to first 
order in the velocity of the W boson and thus not strongly affected by a transverse momentum 
of the W boson induced by gluon radiation at NLO QCD [43] . 

The size of the QCD corrections to the px,; distribution can be reduced by applying a jet 
veto. Table 3 shows the impact of a jet veto on the NLO QCD calculation where we restrict the 
additional parton to a transverse momentum j^xjet (= J3t,w at NLO) < 50 GeV. Of course, the 
jet veto also reduces the size of the photon-induced processes as demonstrated in Table 3. 

In Figures 1 and 2 we show the differential cross sections and the corresponding corrections 
with respect to the transverse momentum px,z and the transverse mass Mx^jy;/, respectively. The 
distributions show the well-known kinks at px,i ~ -Mw/2 and Mx,;/,; ~ -Mw, which are used in 
the W-mass determination. Near these kinks the correction 5qq reaches the order of 10% for 
bare muons and is reduced to about 5% after photon recombination. Near the resonance region, 
the corrections from photon-induced processes are very small. Multi-photon emission, on the 
other hand, reaches the percent level near p^^i « Mw/2 and M^^y^i k, Mw and induces some 
distortion that affects the Mw determination from the shape of the Mt^i^^i distribution. The 
two-loop Sudakov corrections are, of course, completely negligible near the W resonance and are 
thus not displayed in Figures 1 and 2. 
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Figure 1: Lepton-transverse- momentum distribution in LO and corresponding relative 
corrections 5 at the LHC in the SM. 




Figure 2: W-transverse-mass distribution in LO and corresponding relative corrections 
6 at the LHC in the SM. 
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2.6.4 Cross sections and distributions for pp —>■ W+ I'^uiX at the Tevatron 

We also present numerical results for W production at the Tevatron, i.e. a pp collider with a 
CM energy of ^/s = 1.96 TeV. We again use the phase-space cuts of (2.27) and the photon 
recombination procedure specified in the previous section. 

Tables 4/5 and Figures 3/4 display the LO cross section ctq and the various relative corrections 
(5 as a function of the transverse momentum and transverse mass. The electroweak corrections 
are typically of the same size as those discussed for the LHC, and also show the same qualitative 
features. It is obvious that the high-energy Sudakov regime is not phenomenologically relevant for 
W production at the Tevatron, and that the size of the two-loop Sudakov corrections at moderate 
PT,i and Mt,i/iI is no reliable estimate of the theoretical uncertainty from missing electroweak 
two-loop corrections. 

2.6.5 Comparison with existing results on multi-photon final-state radiation 

Multi-photon radiation has also been included in the Monte Carlo programs Winhac [16] and 
Horace [15] within a parton-shower approach in leading logarithmic accuracy. The program 
Horace, in particular, combines multi-photon radiation with the 0{a) electroweak corrections. 
In the following, we will compare our results for multi-photon final-state emission within the 
structure-function approach defined in Section 2.5 with the Horace results as presented in Ref. [1] . 

Adopting the setup of Ref. [1], i.e. the choice of SM input parameters, the identification cuts 
for bare leptons, and the a(0) input-parameter scheme, we find excellent agreement with the 
Horace result for the leading order and the 0{a) corrections, as expected from earlier tuned 
comparisons [12]. 

Figure 5 shows the multi-photon final-state corrections to the distributions in and M'y^uiI 
for the muon final state. Here is the result of our multi-photon correction in the a(0) 

scheme in the setup of Ref. [1]. The band indicates the dependence on the QED scale (2.25), 
with the factor ^ varied between 1/3 and 3. The result of Horace is denoted (^J^ijtf^^- For px,i 
and M'Y^y^l values below and at the Jacobian peak, the relative corrections differ by less than 
0.1%. Above the Jacobian peak the comparison shows larger deviations (about 0.2%) which 
are comparable to the statistical fluctuations of the Horace results. The level of agreement is 
satisfactory, keeping in mind that the structure-function and parton-shower approaches involve 
different approximations. In fact, the difference between these approximations should be viewed 
as a lower bound on the theoretical uncertainty. 

Figure 5 also shows the multi-photon correction (5inuiti_-y, as defined in Section 2.5, which 

in contrast to <5^ui|i_^ contains the formal two-loop effect (a(0) — ac^) ln(m^/(5^) induced by 
changing to a(0) in the leading photonic 0{a) correction in the scheme. Of course, this 
contribution does not reflect a genuine new multi-photon effect but only a different partitioning 
of the total correction into an 0{a) and a higher-order QED part. 

3 Supersymmetric corrections in the MSSM 

Non-standard physics could affect the W-boson cross-section prediction and thus bias the preci- 
sion determination of Standard Model parameters from W-boson observables at hadron colliders. 
To quantify the impact of new physics on the W cross section in a concrete model, we have cal- 
culated the 0{a) electroweak and ©(ctg) strong corrections to pp/pp — W"*" — l^viX within 
the MSSM. 
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pp l+viX at = 1-96 TeV 







^0-00 


75-00 


lOO-oo 


200-00 


300-00 


(T()/j)b 




'^1 7 iQrir'?^ 


0.51112(2) 0.15201 1(6) 0.005G 105(2) 0.00039160(2) 




—2 7(1) 


—5 4(1) 


-7.4(1) 


-9.1(1) 


-14.2(1) 


-18.6(1) 




-1 7(1) 


-2 8(1) 


-4.9(1) 


-6.3(1) 


-10.3(1) 


-13.7(1) 


HI' 


0.020(1) 

V / 


1.5(1) 


2.0(1) 


2.0(1) 


1.4(1) 


0.95(1) 


V /% 

Sudakov' ^ 


-0.0001 


-0.017 


-0.056 


-0.078 


-0.013 


0.20 


^multi-7/% 


n 11 +0-02 


n 34+0-09 


30+° °^ 
'-'•"^'^-0.06 




n 4^+0-09 


0.59^°:}? 




-2.6(1) 


-3.5(1) 


-5.1(1) 


-6.8(1) 


-12.3(1) 


-17.1(1) 




-1.6(1) 


-1.3(1) 


-2.9(1) 


-4.3(1) 


-8.9(1) 


-12.8(1) 


&/% 


11.2(1) 


377(1) 


205(1) 


174(1) 


113(1) 


74.6(1) 


<^QCD / /o 


11.0(1) 


362(1) 


176(1) 


138(1) 


69.9(1) 


34.9(1) 



Table 4: Integrated LO cross sections (Tq for W"*" production at the Tevatron for different 
ranges in pT,/ and corresponding relative corrections 6 in the SM. 



pp l+viX at = 1-96 TeV 



MT,.,//GeV 


50-00 


100-00 


150-00 


200-00 


400-00 


600-00 


do/pb 


706.95(1) 


3.7496(2) 


0.51112(2) 0.152014(6) 0.0056405(2) 


0.00039160(2) 




-2.7(1) 


-5.2(1) 


-6.5(1) 


-8.0(1) 


-12.7(1) 


-16.8(1) 




-1.7(1) 


-3.4(1) 


-4.8(1) 


-6.2(1) 


-10.1(1) 


-13.3(1) 




0.017(1) 


0.028(1) 


0.028(1) 


0.027(1) 


0.018(1) 


0.012(1) 


4udakov/^ 


-0.0001 


-0.017 


-0.056 


-0.078 


-0.013 


0.20 


'^multi-7/% 


11+0-02 




^•■•-0-0.04 


iQ+O-04 


97+0.05 
-0.05 


n 07+0.07 
-0.07 




-2.6(1) 


-4.9(1) 


-6.3(1) 


-7.8(1) 


-12.4(1) 


-16.4(1) 




-1.6(1) 


-3.4(1) 


-4.8(1) 


-6.1(1) 


-10.0(1) 


-13.3(1) 




10.8(1) 


22.1(1) 


19.6(1) 


16.7(1) 


6.2(1) 


-2.8(1) 


j-/i=MT,w /(W 

<^QCD / /o 


10.5(1) 


21.4(1) 


18.5(1) 


15.3(1) 


4.1(1) 


-5.3(1) 



Table 5: Integrated LO cross sections (Tq for W"*" production at the Tevatron for different 
ranges in Mr,^,,/ and corresponding relative corrections S in the SM. 
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Figure 3: Lepton-transverse-momentum distribution in LO and corresponding relative 
corrections 5 at the Tevatron in the SM. 




Figure 4: W-transverse-mass distribution in LO and corresponding relative corrections 
6 at the Tevatron in the SM. 
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Figure 5: Comparison of multi-plioton corrections as obtained in this work with results 
of Horace [1,15] (see text for details). 



3.1 Supersymmetric QCD and electroweak corrections 

The SUSY-QCD corrections comprise gluino-squark loops which contribute to the qq'W vertex 
correction and to the quark wave-function renormalization. 

To obtain the electroweak SUSY corrections, we calculate the complete electroweak 0{a) 
corrections in the MSSM and subtract the SM corrections, so that the MSSM corrections can be 
added to the SM predictions of the previous section without double counting. This procedure 
applies to the vertex, box, and self-energy corrections as well as to the counterterms and Ar. 
The diagrams for the genuine SUSY vertex and box corrections are shown in Figures 6 and 7, 
respectively. The subtraction is only non-trivial for contributions from the Higgs sector where the 
SM and MSSM vertex appears with different couplings. However, for massless external fermions 
the Higgs sector only contributes to the W-boson self-energy and to renormalization constants. 
The mass of the lightest MSSM Higgs boson is used as SM input when we subtract the SM 
contribution to the correction. Hence, for ultimate precision, one would have to calculate the 
SM corrections with the appropriate Higgs mass. 

The diagrammatic calculation is performed in two almost independent ways. Both calcu- 
lations use FeynArts [21] to generate the relevant (MSSM) diagrams. One of the calculations 
then employs FormCalc and LoopTools [22] to perform the algebraic calculation and the loop 
integrals while the other calculation uses a completely independent set of in-house routines for 
both the algebraic calculation and for the numerical evaluation. 

3.2 Numerical results 

3.2.1 Input parameters and setup 

The SM input parameters and the setup of the calculation (input parameter scheme, PDFs, cuts, 
etc.) are chosen as described in Section 2.6.1. 

To study the dependence of the corrections on the SUSY breaking parameters, we show 
results for all the SPS benchmark scenarios [44]. The generic suppression of the genuine SUSY 
corrections turns out to be rather insensitive to a specific scenario. We therefore refrain from 
further restricting the SPS coverage by taking into account recent experimental bounds in favour 
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Figure 6: Classes of diagrams for the additional vertex correetions within the MSSM. 
The neutralinos x^, charginos x, squarks q {q = u, d) and slap tons / represent the different 
possible mass eigenstates. 




Figure 7: Classes of diagrams for the additional box corrections within the MSSM. The 

neutralinos x^, charginos x, squarks q {q = u, d) and sleptons I represent the different 
possible mass eigenstates. 



of a broader scope in the SUSY parameter space. The SPS points are defined by the low- 
energy SUSY breaking parameters which determine the spectrum and the couplings. For the ten 
benchmark scenarios under consideration in this work, this input [45] is tabulated in Appendix A. 

Dependent SUSY parameters, such as Higgs, chargino, neutralino, or sfermion masses, are 
calculated from the SPS input using tree-level relations. Since the impact of the fermion masses 
of the first two generations is negligible, these masses are set to zero in the calculation of the 
corresponding sfermion mass matrices. Following this approach the SUSY corrections do not 
depend on the fermion generations in the partonic process ud — > I'^i'i. In particular, the SUSY 
corrections presented below are valid for both outgoing electrons and muons. 

3.2.2 Corrections to partonic cross sections 

In order to exhibit the typical features of the SUSY corrections it is instructive to first discuss 
the inclusive partonic cross section (Tq evaluated at a fixed partonic CM energy Here, no cuts 
are applied. In Figure 8 we display the LO result and the SUSY-EW and SUSY-QCD corrections 
as a function of the parton energy s/i for the ten different SPS scenarios. 

The SUSY-QCD corrections turn out to be completely negligible for parton energies below 
1 TeV, as the SPS scenarios involve heavy squarks and gluinos which effectively decouple. The 
SUSY-QCD corrections reach the level only when the sum of the gluino and a squark 

mass equals s/S, which typically happens between 1 and 2 TeV for the SPS scenarios. For even 
larger parton energies well above the squark/gluino mass scale the SUSY-QCD corrections turn 
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Figure 8: Total partonic LO cross section (Tq and corresponding relative SUSY-EW and 
SUSY-QCD corrections b as function of the partonic CM energy \fk for the different SPS 
scenarios. 
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negative and become logarithmically enhanced with increasing s. In Figure 8, the corrections are 
only shown up to the phenomenologically relevant region \/I = 2 TeV such that the asymptotic 
behavior is not visible. 

The EW corrections turn on at smaller -s/J and exhibit peaks which correspond to sparticle 
thresholds. Specifically, corrections at the —1% level occur if the sum of a neutralino and a 
chargino mass equals \/I and if the sleptons and squarks in the model are heavy compared to 
this mass scale. The corrections then rise logarithmically with s to positive values. At the mass 
scale of the sleptons and squarks there is an additional negative contribution to the correction 
which is less peaked than the gaugino peaks and which washes out the gaugino peaks if the 
corresponding regions overlap. 

3.2.3 Corrections to hadronic cross sections 

To calculate the hadronic MSSM cross section we apply the same phase-space cuts as for the 
SM analysis. However, since the MSSM corrections are purely virtual, the kinematics of the 
relevant events is as in leading order, i.e. photon recombination is irrelevant and Mx^jy;/ = 2px,z- 
In Table 6, we show the integrated cross section for different ranges in pT,z in analogy to the SM 
analysis in Section 2.6. As expected, the corrections for relatively low pt,/ cuts are negligible at 
the sub-per-mille level. Only in the high-pT,z tail, the EW corrections reach the percent level 
if the SUSY spectrum is light enough. This rise in the corrections in the high-pT,« tail depends 
on the mass scale of the relevant SUSY particles in the loops. The maximum of the corrections 
is reached for the SPS 2 scenario where the gauginos are particularly light and the squarks and 
sleptons are so heavy that their negative contribution becomes effective only at even larger pT,«- 

SPS 2 is also the only scenario for which the EW corrections in the pT,z distribution, as shown 
in Figure 9, almost reach the percent level for px.z < 100 GeV due to the light gauginos. In the 
W-resonance region the corrections are extremely small and flat in pt,/- Only extremely light 
gauginos with masses smaller than Mw could impact the precise determination of the W mass. 
Experimental bounds on the mass of the lightest chargino My.± > 100 GeV [39] rule out such a 
scenario. Note that the SUSY-QCD corrections are multiplied by a factor of 10 in the plot. 

The analogous results for the Tevatron are shown in Table 7 and Figure 10. 

4 Conclusions 

Single- W-boson production is one of the cleanest hadron collider processes and will be used to 
precisely determine Standard Model parameters like the W-boson mass and width. In order 
to match the envisaged high experimental accuracy, it is mandatory to reduce the theoretical 
uncertainty of the cross-section prediction to a level of one percent or better. 

We have studied radiative corrections to single- W-boson production, pp/pp W"*" I'^viX, 
at the LHC and at the Tevatron. We have completed our previous calculation of the 0(oi) elec- 
troweak corrections [9] by including photon-quark scattering processes. These photon-induced 
contributions turn out to be considerable at large lepton transverse momentum. However, they 
do not significantly affect the distribution in the lepton-neutrino transverse mass Mt.i/,/. We 
have furthermore discussed the impact of electroweak effects beyond 0{a), specifically the Su- 
dakov logarithms which arise in the high-energy regime. The leading 0{a^) Sudakov logarithms, 
which we consider as a measure for the electroweak two-loop effects, turn out to be small, below 
5% even for transverse lepton momenta pT,z in the TeV range. Corrections due to multi-photon 
final-state radiation beyond 0{a,) reach the percent level near the W resonance and distort the 
shape of the pt,z and Mr,fc,,z distributions. Given the high experimental accuracy envisaged 
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Table 6: Relative SUSY-EW and SUSY-QCD corrections S in the MSSM for W+ pro- 
duction at the LHC for different ranges in px,;- The corresponding integrated LO cross 
sections ao can be found in Table 1. 
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Figure 9: Lepton-transverse-momentum distribution in LO and corresponding relative 
SUSY-EW and SUSY-QCD corrections 5 for the different SPS scenarios at the LHC. 
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Table 7: Relative SUSY-EW and SUSY-QCD corrections S in the MSSM for W+ produc- 
tion at the Tevatron for different ranges in px,;- The corresponding integrated LO cross 
sections ctq can be found in Table 4. 
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Figure 10: Lepton-transverse-momentum distribution in LO and corresponding relative 
SUSY-EW and SUSY-QCD corrections 5 for the different SPS scenarios at the Tevatron. 
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specifically at the LHC, it is therefore mandatory to control the effects of multi-photon emission 
in the determination of the W mass. 

Finally, to study the impact of new physics on the W cross section in a concrete model, we 
have calculated the 0{a) electroweak and 0{as) strong corrections to pp/pp — W"^ l^^iX 
within the MSSM. The supersymmetric corrections turn out to be negligible in the vicinity of 
the W resonance for viable MSSM scenarios, reaching the percent level only at very high lepton 
transverse momenta and for specific choices of the supersymmetric parameters. Effects from 
virtual SUSY particles would thus not spoil the status of single- W-boson production as one of 
the cleanest Standard Model candles at hadron colliders. 
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A SPS benchmark scenarios 

For the SPS benchmark [44] scenarios discussed in this work we use the low-energy input specified 
in Table 8. The input variables are the ratio of the vacuum expectation values of the Higgs 
bosons giving rise to up- and down-type fermion masses, the mass of the CP-odd Higgs boson, 
M^o, the supersymmetric Higgs mass parameter //, the electroweak gaugino mass parameters 
M\^2-, the gluino mass nig, the trilinear couplings A^-^t.b) the scale at which the DR-input values 
are defined, /LiR(DR), the soft SUSY-breaking parameters in the diagonal entries of the squark 
and slepton mass matrices of the first and second generations M/j (where i = L,R refers to the 
left- and right-handed sfermions, f = q,l to quarks and leptons, and / = u, d, e to up and down 
quarks and electrons, respectively), and the analogous soft SUSY-breaking parameters for the 
third generation Ml^. 
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